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Abstract 

The m-point nonlocal problem for the first order differential equa- 
tion with an operator coefficient in a Banach space X is considered. 
An exponentially convergent algorithm is proposed and justified pro- 
vided that the operator coefficient A is strongly positive and some 
existence and uniqueness conditions are fulfilled. This algorithm is 
based on representations of operator functions by a Dunford-Cauchy 
integral along a hyperbola enveloping the spectrum of A and on the 
proper quadratures involving short sums of resolvents. The efficiency 
of the proposed algorithms is demonstrated by numerical examples. 
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1 Introduction 

The m-point initial (nonlocal) problem for a differential equation with the 
nonlocal condition u(to) + g{t\, . . . ; t p ; u) = uq and a given function g on a 
given point set P = {0 = to < ti < ■ ■ ■ < t p } , is one of the important topics 
in the study of differential equations. Interest in such problems originates 
mainly from some physical problems with a control of the solution at P. For 
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example, when the function g{t\] . . . ; t p ; u) is linear, we will have the periodic 
problem u(t ) = w(ti). Problems with nonlocal conditions arise in the theory 
of physics of plasma [I], nuclear physics [2], mathematical chemistry [3], 
waveguides [1] etc. 

Differential equations with operator coefficients in some Hilbert or Banach 
space can be considered as meta-models for systems of partial or ordinary 
differential equations and are suitable for investigations using tools of the 
functional analysis (see e.g. [HIE]). Nonlocal problems can also be considered 
within this framework [2 [8] . 

In this paper we consider the following nonlocal m-point problem: 

u' t + Au = f(t), te[0,T] 

+ akU ( tk ) = n °> <h <t 2 < ■ ■ ■ <t m <T, 
k=i 

where G R, k = 1, m, f(t) is a given vector-valued function with values in 
a Banach space X, uq G X. The operator A with the domain -D(-A) in a Ba- 
nach space X is assumed to be a densely defined strongly positive (sectorial) 
operator, i.e. its spectrum T,(A) lies in a sector of the right half-plane with 
the vertex at the origin and and the resolvent decays inversely proportional 
to | jar | at the infinity (see estimate ([3]) below). 

Discretization methods for differential equations in Banach and Hilbert 
spaces were intensively studied in the last decade (see e.g. 0, [El [EH [E21 EH 
El US US El HH] and the references therein). Methods from pEl [EH H21 CCS 
El [ISl [El [H] possess an exponential convergence rate, i.e. the error estimate 
in an appropriate norm is of the type 0(e~ Na ), a > with respect to a 
discretization parameter N — y oo. For a given tolerance e such discretizations 
provide optimal or nearly optimal computational complexity [TD]. One of the 
possible ways to obtain exponentially convergent approximations to abstract 
differential equations is based on a representation of the solution through the 
Dunford -Cauchy integral along a parametrized path enveloping the spectrum 
of the operator coefficient. Choosing a proper quadrature for this integral we 
obtain a short sum of resolvents. Since the treatment of such resolvents is 
usually the most time consuming part of any approximation this leads to a 
low-cost naturally parallelizeable algorithms. Parameters of the algorithms 
from [12l [T31 [16] were optimized in [191 EO] to improve the convergence rate. 

The aim of this paper is to construct an exponentially convergent ap- 
proximation to the problem for a differential equation with m-point nonlocal 
condition in abstract setting (pQ). The paper is organized as follows. In 
Section [2] we discuss the existence and uniqueness of the solution as well as 
its representation through input data. An algorithm for the homogeneous 
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problem ([T]) is proposed in section El The main result of this section is theo- 
rem [2] about the exponential convergence rate of the proposed discretization. 
The next section |4] is devoted to the exponentially convergent discretization 
of the inhomogeneous problem. In section [5] we represent some numerical 
examples which confirm theoretical results from the previous sections. 



2 Existence and representation of the solu- 
tion 

Let the operator A in (Tjrj be a densely defined strongly positive (sectorial) 
operator in a Banach space X with the domain D(A), i.e. its spectrum 
lies in the sector (spectral angle) 

£ = {z = Po + re ie : r G [0, oo), \9\ < if < |} . (2) 

Additionally outside the sector and on its boundary i\ the following estimate 
for the resolvent holds true 

The numbers po, if are called the spectral characteristics of A. 
The hyperbola 

ro = { z {£) = Pocosh£ — ib smh.^ : £ G (— oo, oo), b = p tan ip} (4) 

in turn is referred as a spectral hyperbola. It has a vertex at (po, 0) and 
asymptotes which are parallel to the rays of the spectral angle E. 

A convenient representation of operator functions is the one through the 
Dunford-Cauchy integral (see e.g. [51 E]) where the integration path plays an 
important role. We choose the following hyperbola 

r^- = {z(£) = Of cosh£ — ibi sinh£ : £g (—00,00)}, (5) 

as the integration contour which envelopes the spectrum of A. 
For an arbitrary vector uq G D(A m+1 ) the next equality holds 

m+l Ak~l 1 

E + ~ A r lAm+lu « = ( zI - A r^o. (6) 

k=l 

This formula, together with 

A -{m+l) v = J_ I z -(m+l)( zI _ A yl vdz ^ (7) 
27TZ J Tl 
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by setting v = A m+1 u , yields the following representation 



u 



A -(m+l) A m+l UQ = _ 



2ni 



m + 1 Ak-1 



k=l 



-(m+1) 



u dz. 



(zl - A)- 1 A m+1 u d2 



Various useful properties of this representation as well as the next result 
were discussed in 



Theorem 1 Let Uq G D(A m+a ) for some m G N and a G [0,1], then the 
following estimate holds true: 

m+1 A k ~ l 



(zI-A)- 1 -^ 



k=l 



~ \z\ m+1 (1 + \z\) a 11 °"' (9) 

Wa G [0,1], n G D(A m+a ), 

where the constant K depends on a and M only. 

In Ej it was proven that the solution of the problem ([T]) exists and is 
unique provided that one of the following two conditions is fulfilled: 



5^ 1^1 < lj 



or 



i=i 



aAe~ poti < 1 



(10) 



with p from fl2j. 

According to the Hille-Yosida-Phillips theorem [22] the strongly positive 
operator A generates a one parameter semigroup T(t) = e~ tA and the solution 
of ([I]) can be represented by 



u (t) = e - At u(0) + / e- A{t - T) f(r)dr. 



(12) 



Combining the nonlocal condition from (JT]) and ( TT2]) we obtain 



u(ti) =e 



-At,. 



u - y^ftfc^fa 



k=l 



+ 



-A{U-t) 



f(r)dr, 



(13) 



i — 1, m. 



Multiplying these equations by and summing up over i — l,m one gets 



t=i i=i fc=i / 14 x 

/ e- A ^- T )/(r)dr. 



Let us denote 



i=i 



then ( 1T4|) implies 

m m m „^ 

„-_i »_i 



i=l i=l i=l 

m 

Setting 5 = I + oii^~ Ati , we obtain 



i=i 

m „t . 



PL 

= Bu -uo + y^cti / 1 e- A(tl - T) /(r)c/r, 
i=i ^ 



where I is an identity operator. Condition ffTOj) or ffTTj) guarantee the exis- 
tence and boundedness of B~ l (A), therefore we have 



W = u - B- X u + B- 1 V a t / e- A{tl - T) /(r)rfr. 

i=i ^ 

Now, the equations ffT3]) yield the representation for the solution of ([I]) in the 
form u(t) = Uh(t) + Uih(t), where Uh(t) = e~ At B~ 1 u is the solution of the 
homogeneous problem with the initial condition uq and 



u ih (t) = -e- M B- 1 JZa i e - A ^f(r)dr+ e~ A ^f(r)dr (15) 
i=i Jo 



is the solution of the inhomogeneous problem with zero initial condition. 



3 Homogeneous nonlocal problem 

First of all we consider the solution Uhit) = e~ At B~ 1 (A)uo of the homoge- 
neous problem ([1]). Our aim in this section is to construct an exponentially 
convergent method for its approximation. 
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Using the Dunford-Cauchy representation of u^{t) analogously to [2T] we 
obtain 



Uh(t) 



— I e~ zt B- 1 (z)(zI - A)- l u dz 

e -zt 

(zl - A)~ l u dz. 



2 

1 

2txI 



(16) 



1 + Yl a i e ~ 



i=l 



Representation (fT6l) makes sense only when the function e~ zt B~ 1 (z) is ana- 
lytic in the region enveloped by Tj. Let us show, that conditions like flTt)]) or 
(|lip guaranty this analyticity [BJ. 

Actually, the analyticity of e~ zt B~ 1 (z) might only be violated when 

n 

Y aie~ zti = —1, since in this case B~ 1 (z) becomes unbounded. It is easy to 

8=1 

see that for an arbitrary z we have 



-zt k 



k=l 



> 



a; fe • e 



t^Of cosh £ I 



-it^bj sinh £ 



> 



fc=l 
m 



-P0*fc 



fc=l 



> o, 



provided that ffTTj) holds true. 

Using (jHJ) with m = we can modify the representation of Uh(t) as follows: 



u h (t) 



1 
1 

2ttz 



(zl - A)- 1 - -I 



-zt 



Tl l + Y, one- 



(zI-A)" 1 --! 



uodz 
uodz. 



;i7) 



i=l 



After discretization of the integral such modified resolvent provides better 
convergence speed than (TIB"]) in a neighbourhood of t = (see [21] for details). 
Parameterizing the integral ()17|) by © we get 



with 



1 Z" 00 

m*) = ^- y -F(UR, 
j-(u) = f a (uk 
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n 

1 + a^-^ 

i=l 



*(£) \ 



z'(£) = a i sinh £ — i6j cosh £ . 
Supposing wo G -D(A Q ), < a < 1 it was shown in [21] that 

1 



-*(0V 



(^e)i-A)- 1 



Mo | 



<( 



1 + M)K— (—\ e- aitcosh ^- a ^\\A a uo\\, (Gl,i>0. 
ai \aij 



The part responsible for the nonlocal condition in (I18p . can be estimated in 
the following way 



1 4 > a k e 

k=l 



in 



ajfc e 



-a/ cosh(§)i fe 



fc=l 

-1 



Qtfc e 



-pi** 



fe=i 



where a/ > p\. 

Thus, we obtain the following estimate for T{t, £): 

11^,011 < Q(1 + M)K— f — ^ e^* cosh «- a l«l||A a M ||, £ e M, t > 0. 

ai \ai ) 



(19) 



The next step toward a numerical algorithm is an approximation of ( II 8p 
by the efficient quadrature formula. For this purpose we need to estimate 
the width of a strip around the real axis where the function T{t, £) permit 
analytical extension (with respect to £). After changing £ to £ + %v the 
integration hyperbola Tj will be translated into the parametric hyperbola 
set 

r(r/) = {z(u>) = a/ cosh (£ + iz/) — z&j sinh (£ + «/): £ G (— oo, oo)} 
= {z(w) = a(u) cosh£ — ib{v) sinh£ : £ G (— oo, oo)}, 

with 



a(u) = aj cos v + bj sin i/ = W + 6j sin {y + 0/2) 



6(j/) = bi cos i> — ai sin = y a\ + b\ cos + 0/2), 
6/ .0 



cos 



sin 



2 v/^Tfef 2 v /^T6f 
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The analyticity of the function J-"(t, £ + iv), in the strip 



D dl = {(£, v) : £ G (-00, 00), < tZi/2}, 



with some d\ could be violated if the resolvent or the part related to the 
nonlocal condition become unbounded. To avoid this we have to choose d\ in 
a way that for v G (—di/2, di/2) the hyperbola set remains in the right 
half-plane of the complex plane. For v = —di/2 the corresponding hyperbola 
is going through the point (pi,0), for some < p\ < p . For v = di/2 it 
coincides with the spectral hyperbola and therefore for all v G (—di/2, d\/2) 
the set does not intersect the spectral sector. This fact justifies the 
choice the hyperbola T(0) = T/ as the integration path. 

The requirements above imply the following system of equations 



dj cos (di/2) + bj sin (di/2) = p , 

bi cos (d\/2) — ai sin (c?i/2) = b = p tan<p, 

ai cos (— di/2) + bi sin (— d\/2) = p±, 



it lead us to the next system 



a i — Po cos (d\/2) — b sin (di/2) 
bi = Po si n (^i/2) + b cos (di/2), 
2aj cos (di/2) = p + Pi- 



Eliminating ai from the first and the third equations we get 



Po cos di — b sin d\ = p\ 



cos(<ii + ip) 



Pi 



i.e. 



d\ = arccos 




(20) 
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with costp 



Po 



sin <p 



bo 



Thus, for aj, bj we receive 



u 



, a, \ cos arccos 
cos(f + <p) V 
Po — Po 

COS (f 



Tfe /2 + f,/2 



COS <p 



( y + <p 



(21) 



/j, s cos arccos , 
cos (f + <p) \ Vv 7 ^! 
Po L = Po 



pi 



/2 + <p/2 



COS <p 



COS <p 



For aj and 6/ defined as above the vector valued function Tit, w) is analytic 
in the strip with respect to w = £ + %v for any t > 0. Note, that for 
pi = we have d\ = it/ 2 — <p as in [21] . 

Taking into account f!2ip we similarly can write equations for a{v),b{v) 
on the whole interval — 4p < z/ < 4r 

a(z/) = a/ cos z/ + b\ sin z/ = \j p\ + b\ cos + cos(^) 

+ \j Pl + h l sin fy + <pJ sin(z/) = \J pl + bl cos f y + <p - i/Y 

= 6/ cos i/ — a/ sin z/ = \J p^ + b^ sin f y + V 9 j cos(^) 



+ &§ cos + <pj si 



sin z/ 



Po + b sin ( TT + <P - ^ ) > 



Pi < a(z/) < po, &o < &0) < y^o + Po -Pi> 

with d u defined by (j2Dj). 

For the part responsible for the nonlocal condition we have 



t=l 



< 



a,- e 



-*(£,i/)ti 



t=l 



< 



a,-e 



-a(i/) cosh({;)t; 



i=l 



9 



-1 



< 1 - |«i|e- a ^ < 1 - |«i|e-^ < q- 1 = Q, 



i=l / \ i=l 



where z(£, u) = a{v) cosh(£) — ib{y) sinh(£). 
Now, for w G Dd L we get the estimate 



(1 + M)QKJa?(u) sinh 2 £ + b 2 (u) cosh 2 £ 

h(,A -a(u)t cosh e, 

~ V ^ a{u) (a 2 {u) cosh 2 £ + b\v) sinh 2 £)W 2 ) 11 011 

< (1 + M)QK^4 f e- aMtcosh5 - Q|5l |U a n || 
a{v) \a{y)J 

<(l + M)QKt a n(^ + V -u) ( ^ e-^l||^o||, 

Mw G D d . 

Similarly to [23], we introduce the space H p (Dd), 1 < p < oo of all 
vector-valued functions J 7 analytic in the strip 

Da = {z G C : — oo < 9ftz < oo, |Sfcz| < d}, 

equipped by the norm 

= |lim^ (/ aDd(e) 11^)111^1)^ ifl<P<oo, 
\lim e ^o sup zeaDd(e) ||.F(z) || if p = oo, 

where 

D d (e) = {zeC: \Re(z)\ < 1/e, |Im(*)| < d(l - e)} 

and dDd(e) is the boundary of Z?d(e). 

Taking into account that the integrals over the vertical sides of the rect- 
angle -Ddi( e ) vanish as e — > 0, the above estimate for ||.F(t, w)\\ implies 

II^OIlHi^^HA-tiollIC-^a) 

r°° (22) 
+ C + (ip,a)} e- a ^ = C( V ,a)\\A a Uo \\ 1 ; 



with 



C(<p, a) = -[C+(<p, a) + C-(<p, a)}, 
a 

C ± (<p t a) = (l + M)QKtaa^ + <p±^j ( 



2 cosip 



p cos(f + <p±%) 
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Note that the influence of both the smoothness parameter of Uq given by a 
and of the spectral characteristics of the operator A given by tp and po is 
accounted by that fact, that the constant C(<p,a) from fl2Til tends to oo if 
a — > 0, ip — > 7r/2 or p\ — )■ (in this case due to (l20l) d\ — > | — </?). 

We approximate integral (1181) by the following Sinc-quadrature [23| 121]: 

= ^ E F(t,z(kh)), (23) 



k=-N 



with an error 



11^(7*, /i) || = ||ttft(t) 



7 OO 7 

<IMt)-^ E ^(MWJII + II^ E ^MM0)ll 

fc=-oo \k\>N 
I e -7rdi/h 

- 272sinh(7rrf 1 // i )" Jr " Hl(Ddl) 

C(^o;)/i||A Q iio|| A r w , 
+ vy y — 2^ exp[-a 7 tcosh(fc/i) - akh] 

k=N+l 

r \\A a Un\\ f p-T^/h >, 

< J ° A \ — t^-tttt + exp[-ajt cosh ((N + l)/i) - a(N + l)/i] I , 

a { smh (ndi/h) J 

where the constant c does not depend on h, N, t. Equalizing the both expo- 
nentials for t = implies 



h 

or after the transformation 



= a(N + l)h, 



ft = '/^7vTT)- < 24) 



With this step-size the following error estimate holds true 



a 



\Vn(F, h)\\ < -exp (-J^(N + 1) ) \\A a u \\, 



with a constant c independent of t, N. In the case t > the first summand in 
the argument of exp[— a it cosh ((iV + l)/i) — a(iV + from the estimate 
for Hr/A^J 7 , /i)|| contributes mainly to the error order. Setting in this case 
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h = c\\nN/N with some positive constant c\ we remain, asymptotically for 
a fixed t, with an error 



\\r) N (F,h)\\ < C [ e ~-diN/( Cl lnN) +e - Cl a /t 7V/2- Cl aln7Vj \\A a U \\, ( 2 6) 

where c is a positive constant. Thus, we have proved the following result. 

Theorem 2 Let A be a densely defined strongly positive operator and u £ 
D(A a ), a £ (0,1), then the Sine- quadrature represents an approximate 
solution of the homogeneous nonlocal value problem (T7J) (i.e. the case when 
f(t) = 0) and possesses an exponential convergence rate which is uniform 
with respect tot>0 and is of the order (D(e~ c ^) uniformly int > provided 
that h = 0(1/ y/N) (estimate [25])) and of the order O {max{e- ndN/{cilnN \ 
e -c iai tN/2-c ia \nN^ j or mch fi xed t > q provide fatf = Cl \ n N/N (esti- 
mate l[%b)) ). 

4 Inhomogeneous nonlocal problem 



In this section we consider the particular solution ([15]) of inhomogeneous 
problem ([1]), i.e. with f(t) ^ 0. 

Let us rewrite formula f lT5]) in the form 

u ih {t) = u hih (t) + u 2 , ih {t), (27) 

with 

rt m 

ui, ih (t)= e~ A{t - T) f(T)dT, u 2 , ih (t) = -y2otjU2,ihj(t), (28) 

J ° 3=1 

where 



U2,ihAt) = / B-h-^-^fi^dr. (29) 
J o 

We approximate the term Ui^hit) by the algorithm proposed in [21] : 

h N 1 
u 1>ih {t) * Ul , N {t) =— Akh)[(z(kh)I - A)- 1 - ——1} 



where 



2ni z(kh) 

k=~N 
N 

x h VkA^fMt)), 

p=-N 



^k, P {t) = 2 exp ^~2 Z ^ h ^ _tann G°^)]}/ cosn2 (ph)> 
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(30) 



Up (t) = - [1 + tanh (ph)], h = (l/ViV) , 

z(£) = ai cosh£ — ibi sinh£, z'(£) = aj sinh£ — ibi cosh£. 
The next theorem (see [2TJ) characterizes the error of this algorithm. 

Theorem 3 Let A be a densely defined strongly positive operator with spec- 
tral characteristics p , <p and a right hand side fit) G D(A a ), a > for 
t G [0, oo] can be analytically extended into the sector £j = {pe ldl : p G 
[0,oo], < </?} where the estimate 



\A a f(w)\\ < c a e~ SalRewl , weEj 



(31) 



with 5 a G (0, v^Po] holds, then algorithm [3ty) converges to the solution of 
[W with the error estimate 



\£ 



N 



\\v>i,ih(t) - Ui,JV 



(32) 



uniformly in t with positive constants c, c\ depending on a, ip, po an d inde- 
pendent of N. 

Let us construct an exponentially convergent approximation to the term 
U2,ih- Using the representation of the operator functions by means of the 
Dunford-Cauchy integral for the j-th summand of Ui^h we get 

tWO = r f z- z{t+h - s) B~\z)[{zI - AY 1 - -I}f{s)dzds 

J0 ATM JTt z 



2ixi 



mi -a)- 1 



(33) 



f(s)dsz'(OdS, 



z(£) =ai cosh£ — ibi sinh£. 

Replacing here the first integral by the Sinc-quadrature with h = O (iV -1 / 2 ) 
we obtain 



U2,ih,j(t) ~ U 2 ,N(t) 



h 



2-ki 



X 



N 

£< 

k=-N 



-z(kh)t i 



s!{kh)B- 1 {z{kh)) 



(34) 



(z(kh)I - A)^ 1 



z{kh) 



f 



where 



fk,j 



-z(kh)(tj—s) 



i- s, f(s)ds, k = -N,...,N,j = l,n. 



(35) 
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In order to construct an exponentially convergent quadrature for these 
integrals we change the variables by 

s = |(l+tanh£), (36) 

and obtain the improper integrals 

/oo 
MhiOdt, (37) 
-oo 

with 

= Q \ 2 exp[-z(kh)t,(l - tanh0/2]/(tj(l + tanhfl/2). (38) 
2 cosh £ 

Note that equation (I3"6"j) represents the conformal mapping w = ip{z) = 
tj[l + tanhzJ/2, z = (f>(w) = |ln^jp of the strip D v onto the eye-shaped 
domain domain A v [23j . On the real axes the integrand can be estimated by 



|.Ffc(tj,OII < -Tl - exp[-o f cosh(A;/t)tj(l - tanhg)/2] 
2 cosh 4 

x 11/(^.(1 +tanh0/2)|| 



(39) 



<2t ie - 2 '€l 11/(^(1 + tanli0/2)l|. 
The following result from [21] is needed to justify a Sinc-quadrature for 



Lemma 4 Let the right-hand side f(t) in (TJP for t £ [0, oo] can be analyti- 
cally extended into the sector Hf = {pe 101 : p £ [0, oo], \9i\ < tp} and for all 
complex w £ S/ we have 

\\f(w)\\<ce~ s \ Rew \ (40) 

with 5 £ (0, y/2po], then the integrand J-k(tj,£) can be analytically extended 
into a strip Dd i: < d\ < <p/2 and belongs to the class H 1 (Drf 1 ) with respect 
to where po? V 9 are ^ e spectral characteristics of A. 

Under the assumptions of Lemma H] we can construct the following expo- 
nentially convergent Sinc-approximation to fkj : 

N 

fkj ~ fk,j,N = h ^ ^k,p,jf(^p,j), (41) 
p=-N 
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where 



f^k,p, 



P.J 



tj exp{ — -^z(kh)[l — tanh (ph)]} 
2 cosh 2 (ph) 

^[l+tanh(p/i)], h = 0(l/VN), 



z(£) = (ij cosh£ — ibj sinh£. 

Substituting fT4Tl) into fl34"|) we get the following algorithm to compute an 
approximation u 2 j^n to u 2 ,ih,j 



h 



N 



u 2 , ihd (t) & u 2d>N (t) = — E ^ z{kh)t ^{kh)B-\z(kh)) 



k=-N 



[z(kh)I-A)' 1 



1 



(42) 



p=-N 



z(kh) 

We represent the error in the form 

£ N (t) = u 2 ,ih,j(t) - u 2)j , N (t) = r 1)N (t) + r 2)N (t) 

where 

ri,N(t) = u 2 , ihJ (t) - u 2>N (t), 

?"2,7v(t) = U 2)N (t) - U 2>jtN (t). 

Using estimate (1251) (see also Theorem [2]) we obtain the estimate for r x >A r(t) 



(43) 
(44) 



\\n,N 



N 



2ni 



F A(t + tj - s, kh) \ f(s)ds 



k=-N 



(45) 



< — exp 

a 



nda 



(N + l) 



\A a f(s)\\ds, 



where F^(t, £) is the operator defined in sectional Due to (JH]) for m = we 
have for the error r 2 n 



A y e -^ kh ) t z '(kh)B(z(kh)) 

1 



N 



k=-N 



X 



< 



(z(kh)I -A)- 1 
h(l + M)QK 



R 



z(kh) 

-^Wz'ikh)] 



(46) 



2tt 



E 

k=-N 



\z(kh)\ 



1+a 



\A a R 



k,j\ 
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where 



The estimate fl39l) yields 



Rk,j fk,j fk,j. 



N ■ 



11^(^,011 < 2t,e- 2 ^||^/(^(l +tanh0)||. 



(47) 



Due to Lemma H] the assumption ||A Q /(w)|| < c a e Sa \ Rcw \ \Jw G Ej guaran- 
tees that A a f(w) e B}(D dl ) and A a J T k (t j ,w) G H^AfJ- These two condi- 
tion turns us to the situation as in proposition of Theorem 3.2.1, p. 144 from 
[23] with A a f(w) instead of / which implies 



\\A a R k>j \\ = \\A a (f k;j -f k;j>N ) 



/oo °° 
A a 7 k {t 3) i)di-h A a T k {t 3 M) 
-co „ 



< 



k=— oo 
e ndx/h 



h A a T k (tj,kh) 



\k\>N 



+/i ^ 2t j e- 2|fch| ||A a /(^-(l + tanhA;/ i ))|| 

|fc|>iV 

<ce- 2 ^/ ft ||^/(t„t,)|| H1 (D dl ) 
+/i ^ 2t ie - 2|fch| c Q exp +tanhfc/i)| 



|fc|>iV 



i.e., we have 



\\A a R k J < ce~ c ^, 



(48) 



where positive constants c a , 5 a , c, C\ do not depend on t, N, k. Now, 
takes the form 



N 



\Mt)\\ 



2ni 



x 



e- z{kh)t z'(kh)B(z(kh)) 
1 



k=-N 

(z(kh)I-A)- 1 



z{kh) 



Rk,j 



(49) 



with Spf(t) = ^2 k= _ N h |c \z( kh }\\+T n ■ Using the estimate (4.8) from [21] and 

\z(kh)\ = J a] cosh 2 (kh) + b] sinh 2 (fc/i) 
> a 7 cosh(A;/i) > a 7 e |feh| /2, 



< ce- ClVW S N (t), 



\e- z( - kh ^z'(kh)\ 



(50) 
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the last sum can be estimated by 
r - r- 

\S N (t)\ < —= e - Q|fc/VF| < c / e~ at dt < c/a V t G [0, oo). (51) 
VN k <f N J-VN 

Taking into account (|48|) . ( ]5"Tj) we deduce from (l49i) 

||r 2lJV (*)|| <ce~ c ^. (52) 
The next assertion now follows from (143)) . (1451) and (I52p . 

Theorem 5 Lei £/ie assumptions of theorem [3] /ioW. T/ien algorithm fjjfy 
converges uniformly with respect to t and moreover the following error esti- 
mate holds true: 

||M*)II = IKi^(t) - u 2d>N {t)\\ < ce~ ClVW , (53) 
with positive constants c, c\ depend on a, ip, po and independent of N. 
Now we can use approximation (I30p for the every summand and get: 

u 2 ,ih(t) = / B- l e-^ t+u -^f{r)dr « J2 u ^At) = u 2 , N (t)- (54) 

i=l 7 ° j=l 

Theorem [5] guarantees that the error this approach will be bounded by: 



\u 2 , ih (t) ~ U2, N (t)\\ < J2 ce ~ C1 ^ ^ c 2 e" C1 ^. (55) 



Thus, the approximations ( I23|) together with (l30|) and (l54|) represent an 
exponentially convergent algorithm for the problem ([1]). 

5 Numerical examples 

Example 6 We consider the homogeneous problem (QJ) with the operator A 
defined by 

D(A) = {u(x) G H 2 (0, 1) : u(0) = u(l) = 0}, 

Au = -u"(x) \/u G D(A). [ ' 

The initial nonlocal condition reads as follows: 

u(x, 0) + 0.5u(x, 0.2) + 0.3u(x, 0.4) = (1 + 0.5e~^ - 2 + 0.3e _7r2 °- 4 ) sin(Trx), 
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with u = (1 + O^e-^ - 2 + 0.3e~ n2oA ) sm(vrx) G ITie ezaci solution of 

the problem is u(x,t) = e~ n *sin(7ra;). It is easy to find that 

z + sin(7ra) = _A_i. 

Calculations with the algorithm [JW above has been performed in Maple 
with h = N' 1 ' 2 . The error at x = 0.5, t = 0.3 is presented in Table U\ and 
clearly exhibits an exponential decay according to the theoretical estimates. 



N 


En 


4 


.29857983847712589e-l 


8 


.41823888073604986e-2 


16 


.11258594468208641e-2 


32 


.10042178166563831e-3 


64 


.28007158539828452e-5 


128 


.2098826601399176e-7 


256 


.1858929920173152e-10 


512 


.856837124351510e-15 



Table 1: The error for x = 0.5, t = 0.3. 

Due to Theorem® the error should not be greater then En = O (e~ c ^^ . 
The constant c in the exponent can be estimated using the following a-posteriori 
relation: 

c = hx(^-){y/2- lr'N- 1 ' 2 = In (fi N ) (y/2 - l)" 1 ^ 1 / 2 . 

The numerical results are presented in the Table® for this estimation show 
that the constant can be estimated as c ~ 1.5 when N — > oo. 

The next example deals again with a homogeneous problem but in this 
more realistic case the resolvent of A on the element Uq cannot be calculated 
analytically. 

Example 7 We consider the homogeneous problem (Tjp with the operator A 
defined as in l[5b}) and with the following initial nonlocal condition: 

u(x, 0) + u(x, 0.5) = x ln(x), 



18 



N 



16 

32 

64 

128 

256 



2.372652515388745588587496 
1.120 1487327954495 15627946 
1.458741976765153165445005 
1.527648924601130131250452 
1.476794596387591759032900 
1.499935011373075736075927 
1.506597339081609844717370 



Table 2: The estimate of c 



where Uq = x\n(x) G A a , a < 1/2. In this case the resolvent can be repre- 
sented using the Green function 

(zl — A) 1 u = I z + — r I x\n(x) = / G(x, s)s\n(s)ds, 
V dx 2 J J 



G(x, s) 



sm(xy/z) sin((l — s)y/z) x < s, 
y/zsm(y/z) 1 sm(sy/z) sin((l - x)y/z) x>s 



where the integrals were computed by exponentially convergent Sine- quadrature 
(see e.g. JEEj) using Maple. The results for x = 0.5, t = 0.3 are presented in 
Tabled 



N 


u(x, t) 


4 


-.241535790017043e-l 


8 


-.228401191029108e-l 


16 


-.194273285627507e-l 


32 


-.192905848633180e-l 


64 


-.192911920318628e-l 


128 


-.192907849909929e-l 


256 


-.192907820740651e-l 



Table 3: Values of the solution u(x,t) for x = 0.5, t = 0.3. 

It can be easily seen that the number of stabilized digits increases according 
to the theoretical prediction by Theorem 

Example 8 Let us consider the inhomogeneous problem (TjJ) with the same 
A defined by (f56j). and the nonlocal condition 

u(x, 0) + 0.5u(x, 0.2) = (1 + 0.5e a2 ) sin(vrx), 
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For f(t,x) at the right-hand side of the equation (QJ) we set 

f(x,t) = (1 + 7r 2 )e'sin(7rx). 

The exact solution of the problem is u(x,t) = e*sin(7nr). We have used the 
algorithm defined by (d2P, $HW> ( f^jP an d implemented in Maple. For x = 
0.5, t = 0.3 the results presented in Table\4\ and are again in good agreement 
with the theoretical predictions. 



N 


En 


4 


.202211483120243 


8 


.726677678737409e-l 


16 


.138993889900620e-l 


32 


.143037059411419e-2 


64 


.554542099757830e-4 


128 


.532640823981411e-6 


256 


.730569324317506e-9 


512 


.648376079810788e-13 



Table 4: The error for x = 0.5, t = 0.3. 
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